1/1 


UNCLASSIFIED 


ELECTROMAGNETIC  DESCRIPTION  OF  THE  THO  STREAM 
INSTABILITV  IN  A  RELATIVISTIC  ELECTRON  BEAM(U>  NAVAL 
SURFACE  UEPPONS  CENTER  SILVER  SPRING  MD  H  S  UHM  MAV  84 
NSHC/MP-84-236  F/G  28/8  NL 


AD- A 179  864 


.»  J*  V  ~-P  "J «  ”  V  -  TT" 


rvn- *  v«  ITT^.’V ^ "WV,  l" r r r  r*  »V»V ^T: vr; y_,y.,,vry."Vl  V 


ft 


IfflLflLLCOPJf 

ELECTROMAGNETIC  DESCRIPTION  OF  THE  TWO 

STREAM  INSTABILITY  IN  A  RELATIVISTIC 
ELECTRON  BEAM 


I  r 

ft.. 


HAN  S.  UHM 


NSWC  MP  84-236 

PLASMA  PHYSICS  PUBLICATION  NO.  84-4 
MAY  1984 


NAVAL  SURFACE  WEAPONS  CENTER 

White  Oak,  Silver  Spring,  Maryland  20910 


V  ^ 
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INSTABILITY  IN  A  RELATIVISTIC  ELECTRON  BEAM 
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The  influence  of  electromagnetic  effects  on  the  two  stream  instability  in 
a  relativistic  electron  beam  propagating  through  a  collisionless  plasma 
channel  is  examined,  and  closed  algebraic  dispersion  relation  for  the  complex 
eigenfrequency  o>  is  obtained.  It  is  shown  that  electromagnetic  effects  can 
have  a  strong  stabilizing  influence  on  the  two  stream  instability  in  a 
relativistic  electron  beam,  drastically  enhancing  the  critical  beam  current 

for  instability. 
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One  of  the  most  basic  instabilities  that  characterize  a  relativistic 
electron  beam  propagating  through  a  collisionless  plasma  is  the  two  stream 
instability,1"3  which  resulted  from  the  relative  drift  motion  between  the  oeam 
electrons  and  the  background  plasma  particles.  Although  the  two  stream 
instability  is  relatively  familiar  in  the  plasma  physics  community,  most  of 
the  previous  studies1’2  on  this  instability  have  been  limited  to  one 
dimensional  calculations.  In  recent  literature,3  Bogdankevich  and  Rukhadze 
investigated  the  two  stream  instability  of  a  relativistic  electron  beam, 
including  the  finite  radial  geometry  effects  on  stability  behavior.  Thus, 
they  were  able  to  determine  the  limiting  beam  current  due  to  the  two  stream 
instability.  However,  their  calculation  has  been  based  on  the  electrostatic 
approximation.  Although  this  is  a  reasonable  approximation  for  a  mildly 
relativistic  electron  beam,  a  significant  modification  to  the  stability 
behavior  is  expected  for  a  ultrarelativistic  beam;  where  the  electromagnetic 
effect  often  plays  an  important  role.  In  this  letter,  the  influence  of 
electromagnetic  effects  on  the  two  stream  instability  is  investigated  for  a 
relativistic  electron  beam  propagating  through  a  collisionless  plasma  channel. 

The  analysis  in  this  paper  is  carried  out  within  the  framework  of  a 
macroscopic  cold  fluid  model  assuming  that  the  beam-plasma  fluids  are  immersed 
in  a  uniform  axial  magnetic  field.  For  purposes  of  analytic  simplification, 
the  stability  analysis  is  specialized  to  the  case  of  a  sharp-boundary 
equilibrium  in  which  the  beam  and  plasma  channel  have  rectangular  density 
profiles,  i.e.. 
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where  the  subscript  j  =  b,  i,  e  denote  beam  electrons,  plasma  ions  and 
electrons,  respectively,  r  =  Rc  is  the  radial  location  of  a  grounded 
conducting  wall.  Since  the  unstable  mechanism  of  the  two  stream  instability 
is  mostly  due  to  the  fluctuations  of  the  axial  electric  field,  the  stability 
analysis  in  this  article  is  restricted  to  perturbations  of  the  transverse 
magnetic  mode  polarization.  In  addition,  it  is  further  assumed  that  the 
perturbations  are  axisymmetric  (3/30  =  0).  This  is  a  reasonable  assumption 
because  the  axisymmetric  perturbation  is  the  most  unstable  mode  for  this  kind 
of  instability.  Combining  all  of  these  assumptions  and  restrictions,  it  is 
straightforward  to  show  that  the  eigenvalue  equation  of  the  two  stream 
instability  is  expressed  as 
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where  6E  (r)  is  the  perturbed  axial  electric  field,  5p(r)  =  l  e.6n.  and 
^  <«  J  J 

a  n  r\  a  J 

<5J  (r)  =  l  eilrj  <5n.  +  nVsv-  )  are  the  perturbed  charge  and  the  axial 
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J  0 
component  of  the  current  densities,  respectively,  e.  is  the  charge,  n.  is  the 

J  J 

equilibrium  density,  v9  is  the  mean  equilibrium  axial  velocity,  6n.  is  the 

J  fc  J 
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perturbed  density  and  6V.  is  the  perturbed  axial  velocity,  of  species  j,  w  is 

J  ^ 

the  complex  eigenfrequency  and  k  is  the  axial  wavenumber. 


Consistant  with  polarizations  of  the  two  stream  instability,  the 
transverse  component  of  the  perturbed  electromagnetic  field  is  neglected  on 
the  present  analysis.  This  approximation  tremendously  simplifies  the 
linearized  fluid  calculation  of  the  equation  of  motion  and  continuity  equation 
for  the  beam  and  plasma  fluid  element.  Thus,  after  a  straightforward 
calculation,  the  perturbed  axial  velocity  is  expressed  as 


iXfo 


■sr.'.aa 


vv> 


6V 


JZ 


J 

~T~ 

Tjmj 


5E 


kV 


z* 


JZ 


(3) 


where  iiij  is  the  rest  mass  of  particles  of  species  j.  Similarly,  the  perturDed 
density  is  given  by 
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More  detail  calculation  including  the  transverse  component  of  the 
electromagnetic  fields  are  currently  under  investigation  by  the  author  and 
will  be  presented  elsewhere.  Note  from  Eqs.  (3)  and  (4)  that  both  the 
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perturbed  density  6n.  and  axial  velocity  6V-  profiles  can  eventually  be 
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expressed  in  terms  of  the  perturbed  axial  electric  field  6EZ.  Substituting 
Eqs.  ( 3 }  and  (4)  into  Eq.  (2)  and  carrying  out  a  straightforward  calculation, 
the  eigenvalue  equation  in  Eq.  (2)  is  expressed  as 
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where  ur.(r)  =  4nnV(r)e^/m .  is  the  nonrelativistic  plasma  frequency-squared  of 

r  J  J  J  J 

2  *2 

particles  of  species  j.  Note  that  urjr)  =  w  .  =  const.,  for  0  <  r  <  Rb  and 


zero,  otherwise. 


For  the  present  purposes,  the  eigenvalue  equation  (5)  is  solved  for  long 
wavelength  perturbations  characterized  by 
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Within  the  context  of  Eq.  (6),  the  physically  acceptable  solution  to  Eq.  (5) 
can  be  approximately  given  by 


1,  0  <  r  <  R.  , 

6E(r)  -a  D  (7) 

'  -tn{r/Rc)/tn(Rc/Rb) ,  Rb  <  r  <  R£ , 

where  a  is  a  constant  and  Rc  >  >  Rb  is  assumed  in  order  to  satisfy  Eq.  (6) 
self-consistently.  Multiplying  Eq.  (5)  by  r6Ez(r)  and  integrating  over  r  from 
r  =  0  to  r  =  Rc,  it  is  straightforward  to  obtain  the  dispersion  relation 
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where  use  has  been  made  of  Eq.  (7).  In  Eq.  (8),  or  =  4ne  n^/m  is  the 
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nonrelativistic  plasma  frequency-squared  for  beam  electrons,  fg  =  ne/nb.  and 
the  ion  contributions  are  neglected.  In  obtaining  Eq.  (8),  it  is  assumed  that 
the  beam  electrons  have  the  same  axial  velocity  Bbc  which  is  related  to  the 

o  p  "1 

mass  ratio  yb  by  =  (1  -  Bb)  .  Axial  velocity  of  the  plasma  electrons  is 
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assumed  to  be  zero.  Within  the  context  of  Eq.  (6),  the  term  k  -  ur/c  in  the 
right-hand  side  of  Eq.  (8)  can  be  neglected.  However,  I  keep  this  term  for 
the  clarity  in  the  subsequent  stability  discussion,  without  affecting  the 
final  results  in  this  letter. 

Two  points  are  noteworthy  from  Eq.  (8).  First,  it  is  obvious  that  the 
term  2/Rb*n( Rc/Rb)  in  Eq.  (8)  represents  the  influence  of  finite  radial 
geometry.  Somehow,  neglecting  this  term  recovers  the  familiar  one  dimensional 
dispersion  relation  for  the  two  stream  instability.  Second,  the  influence  of 
electromagnetic  effects  are  incorporated  into  the  terms  proportional 
to  u  /c  in  Eq.  (8).  Neglecting  these  terms  gives  the  dispersion  relation 
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for  the  electrostatic  approximation.  Paralleling  the  standard  method 


developed  in  the  previous  literature,  ’  the  critical  beam  current  for 


i  istability  is  analytically  found  from  Eq.  (9)  and  is  given  by 
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for  the  electrostatic  approximation,  where  IA  *  mcJ/e  -  17000  ampere  is  the 


Alfven  critical  current.  The  electron  beam  with  current  below  the  critical 


current  is  stable  whereas  the  beam  with  current  above  the  critical  current  is 


unstable.  Note  from  Eq.  (10)  that  the  critical  current  is  inversely 


proportional  to  the  plasma  electron  density(f„). 


In  order  to  demonstrate  the  influence  of  electromagnetic  effects  on 


stability  behavior,  Eq.  (8)  is  numerically  solved  for  complex  eigenf requency 
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u  =  o>r  +  iui.j,  determining  the  stability  boundaries  in  (w ,yb)  parameter  space 


for  specified  values  of  fe.  The  critical  current  Icm  of  the  electromagnetic 


calculation  is  obtained  from  these  stability  boundary  curves.  For  convenience 


in  the  subsequent  analysis,  we  define  the  electromagnetic  current  enhancement  £ 
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which  is  a  function  of  yb  and  fe  in  general.  Shown  in  Fig.  1  is  the 
electromagnetic  current  enhancement  £  versus  the  beam  energy  yb  obtained  from 
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Eqs.  (8)  and  (11)  for  Rc/Rb  =  7.4,  and  fg  =  1  and  fg  =  10.  Obviously  from 
Fig.  1,  the  critical  current  for  the  electrostatic  approximation  is  valid  only 
for  a  mildly  relativistic  electron  beam.  The  electromagnetic  current 
enhancement  £  increases  drastically  with  increasing  values  of  energy  Yb  and 
plasma  density  fg,  thereby  exhibiting  a  substantial  enhancement  of  the 
critical  current  by  the  electromagnetic  effects  for  a  relativistic  beam. 


It  is  instructive  to  rewrite  Eq.  (8)  by 
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which  is  plotted  in  Fiy.  2.  When  the  value  of  the  local  minimum  point  labeled 

P  is  less  than  unity,  Eq.  (12)  has  six  real  roots  and  the  system  is  stable. 

It  is  noted  from  Eq.  (12)  that  when  $b  approaches  to  unity,  the  strength  of 
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the  destabilizing  term  proportional  to  u>p/Yb  reduces  drastically  in  comparison 
with  the  stabilizing  term  proportional  to  2c£/Rbtn(Rc/Rb) .  In  Fig.  2,  the 
influence  of  vertical  line  w  =  kc  drastically  pull  down  the  local  minimum 
point  P  as  the  line  w  =  kebc  approaches  to  u  =  kc,  eventually  providing  strong 
stabilizing  effects.  The  stabilizing  influence  of  electromagnetic  effects 

A 

originates  from  the  perturbed  axial  current  density  60b(r)  in  Eq.  (2). 
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